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ABSTRACT 



A detailed study of the blue supergiant UIT005 (B2-2.5Ia+) in M 33 is presented. 
The results of our quantitative spectral analysis indicate that the star is a very luminous, 
logL/L0 ~ 5.9dex, and massive, M~5OM0, object, showing a very high nitrogen- 
to-oxygen ratio in its surface (N/0~8, by mass). Based on the derived Mg and Si 
abundances, we argue that this high N/0 ratio cannot be the result of an initial low 
O content due to its location on the disk of M33, a galaxy known to present a steep 
metallicity gradient. In combination with the He abundance, the most plausible inter- 
pretation is that UIT 005 is in an advanced stage of evolution, showing in its surface N 
enrichment and O depletion resulting from mixing with CNO processed material from 
the stellar interior. A comparison with the predictions of current stellar evolutionary 
models indicates that there are significant discrepancies, in particular with regard to 
the degree of chemical processing, with the models predicting a much lower degree 
of O depletion than observed. At the same time, the mass-loss rate derived in our 
analysis is an order of magnitude lower than the values considered in the evolutionary 
calculations. Based on a study of the surrounding stellar population and the nearby 
cluster, NGC588, using HSTAVFPC2 photometry, we suggest that UIT 005 could be 
in fact a runaway star from this cluster. Regardless of its origin, the derived param- 
eters place the star in a region of the Hertzsprung-Russell diagram where Luminous 
Blue Variables are usually found, but we find no evidence supporting photometric or 
spectroscopic variability, except for small Ha changes, otherwise observed in Galactic 
B-type supergiants. Whether UIT 005 is an LBV in a dormant state or a regular blue 
supergiant could not be discerned in this study. Subsequent monitoring would help us 
to improve our knowledge of the more massive stars, bridging the gap between regular 
and more exotic blue supergiants. 
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1. Introduction 



There is no doubt that massive stars, those bom with masses above ^SM©, played, and 
still play, a fundamental role in shaping the universe as we observe it. Unmistakable signatures 
of their prese nce are imprinted i n the spectra of star forming galaxies, both nearby and at 



high redshift (|Steidel et al 



19961). It is in fact because of their very existence that the physical 



properties of these distant gal axies can 



?e inv estigated, through directly related stellar indicators 



such as stellar wind features (|Rix et al 



20041) . or indirectly related stellar signatures, such as 



the io nizing radiation fie 



lines ( Kobulnickv & Koo 



d reprocessed by the gas and radiated in a few forbidden emission 



200C 



Pettini et al 



200 ih . Throughout their entire life, from the Main 



Sequence to their end in a supernova event, massive stars drive, to a significant extent, the 
evolution of galaxies via radiative, kinetic and chemical feedback. Understanding their physical 
properties and their evolution in different environments is therefore a relevant issue for many 
different fields in modem Astrophysics. 



A proper understanding of these objects needs to be assessed by testing the level of agreement 
between the different theories involved, namely the atmosphere and interior structures. Predictions 



under different conditions have to be compare d at any given evo 



been a trivial task for the last two decades (see 



Maeder & Meynet 



utionary stage. This h as not 



2000l or the review by 



Herrero 



20081), and although in the last years an increasingly better agreement has been reached, there 



are still some pending issues, such as the efficie ncy of rotational rnixing or the presence of 



other phenomena a ffecting the stellar evolution (|Tmndle & Lennon 



2005 



Hunter et al. 



2008 



Maeder et al. 



2009b 



Based on observations made with the William Herschel Telescope operated on the island of 
La Palma by the Isaac Newton Group in the Spanish Observatorio del Roque de los Muchachos of 
the Instituto de Astroffsica de Canarias. 
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Particularly important in this context is the role that stars in nearby galaxies can play, where 
very different conditions can be found. Only a relatively small number of Galactic massive 
stars are available for these comparative studies, since the massive stellar population is basically 
concentrated in the Galactic plane, where the thick clouds of dust keep them hidden. This small 
number makes the observation of short-lived evolutionar y phases, whic h could be crucial as they 
are affected by strong episodes of mass-loss and mixing (lMasseyll2003b . even more difficult. The 
problem is less severe in nearby galaxies, such as M 33, where we can have access to a much 
larger number of massive stars and objects in these short-lived phases. 



We came across a very interesting object in our spectroscopic survey of massive blue 
supergiants in M33. It seemed quite pr obable that UIT005 (RA: 0l''32'"42.89^ Dec: 



-h30°38'47.2", J2000, 



Massev et al. 



19961) could be in an advanced phase of evolution, one of 
those very short-lived phases and close to the LBV stage. The aim of this paper is to test the level 
of agreement between the results obtained from the model atmospheres and the theory of stellar 
evolution in such object. This will provide us with clues about the consistency between these 
models in an environment and evolutionary stage different from the most studied cases. 

The rest of this paper is organized as follows; we describe in f|2]the observations and data 
reduction steps. The different aspects of the spectral analysis are presented in ^and the results 
are discussed in ^ Finally, our conclusions are presented in ^ Throughout this work, metal 
abundances are expressed as A(X) = \og{N{X) / N(}i)) + 12. Regardi ng the host galaxy, we 



adopt a distance modulus of /i = 24.56 ±0.1 mag (|Freedman et al 



2001 



), and R9S =28.77^ a 



position angle of the line of nodes ^ = 22° and inclination / = 54° (IVila-Costas & Edmunds 
with the galactic center located at RA0l''33'«50.9" Dec -i-30°39'36" (J2000). 



19921) . 
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2. Observations 



2.1. Photometry 

Photometric information for UIT 005 is available in several bands from diverse sources 
Its designation as UIT 005 originates from the work by 



Massev et al, 



(|1996|) , where, along with 



Massev et al, 



( 2006 ) 



fluxes in two ultraviolet bands, U-, B- and V-band photometry was provided 
published new U-, B-, V-, R- and I-band photometric measurements in a recent epoch. In this 
work, UIT 005 can be identified as the source J013242. 9+303847.0. A much fainter source, 
J013242. 86+303847.0, is almost coincident with our star (see below). 

The nearby giant H II region NGC 588 and its ionizing cluster were targeted with the 
Wide-Field and Planetary Camera 2 (WFPC2) on board of the Hubble Space Telescope (HST) 
during Cycle 4 (program 5384), using a suite of 5 different filters: F170W, F336W, F439W, 
F469W and F547M. Given its close proximity (both objects are separated by ~30", see Fig.fl]), 
UIT 005 was imaged by one of the three detectors of the Wide-Field Camera. The fainter source 
mentioned before is spatially resolved in the H ST images. We have carried out the photometry 



of these data using HSTphot (|Dolphin 
transfo rme( 
09951) and 



2000ar) . Magnitudes in the WFPC2 filter system can be 



transfo rme d to the Iqhnson- Cousins system by applying the recipes presented by 
3 



Holtzman et al. 



Dolphinl(l2000bh 



The contamination of 



Massev et a 



estimated, either by using iMassey et al. 



(119961) photometry due to the faint source can be 



(|2006l) or our WFPC2 photometry. Both sets indicate 
a contribution to the combined flux of ~2% for the U, B and V bands. This would not the 



Massev et al. 



(l2006h . the faint 



case for the R-band where, according to the values published by 
source is only ^0.2 mag fainter than UIT 005, which translates into a flux ratio of 1:0.78. 
This is unexpected, given the almost constant flux ratio in the bluer bands, and the fact that 
UIT 005 was detected in the I-band, whilst the faint source was not. Moreover, there is no 
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indication of a second source in the red spectra (see later); if the R-band magnitudes are right, 
then we would overestimate the continuum by almost a factor 2, which would severely affect 
the measured equivalent widths of the metal lines around Ha. We conclude that the published 
R-band magnitude of the faint source is affected by some unknown problem, and that it is not 
representative of its real value. 

If we disregard for a moment the possibility of photometric variability for UIT005, and 
combine the WFPC2/HST photometry with the previous ground based values, we find for 
the mean values V= 16.79±0.05, B-V = 0.014±0.014 and U-B = -0.942 ±0.099 mag. The 
uncertainties associated with these values are an order of magnitude larger than the individual 1 (7 
uncertainties quoted in the different photometric measurements. The small scatter in the V-band 
magnitude and almost constant B-V color suggest a low level of photometric variability, within 
the expectations for blue supergiant stars. At the same time, the scatter in U-B is a little bit higher, 
indicating perhaps photometric variations in the U-band. 

The combination of the different photometric observations could be affected in principle by 
the time elapsed between the measurements. Most unfortunately, the possibility of time variability 
could not be investigated with any of the recent surveys of M 33; UIT 005 is located in a region of 



the galaxy that was not monitored either by the DIREC T survev (iMacri et al. 



more recent WISE observatory M33 variability proje ct (IShporer & Mazeh 



to the deep survey of the whole galaxy introduced by 



Hartman et al 



20011) nor by the 



20061). With regard 



(120061) . UIT 005 appears 



saturated in the images (D. Bersier, private communication). 



Regarding 



from 



Cutri et al. 



onger w avelengths, UIT 005 can be matched to the object 2MASS 01324287-1-3038468 
(120031) ; the star was detected only in the J-band. Full c overage in near-IR bands 



Cioni et al 



(|2009h ; the agreement in the 



have been obtained from the IR survey presented by 
J-band magnitudes of these two surveys is very good. Our object is also cross-identified as 
the Spitzer Space Telescope source SSTM3307 J01324289-I-303847.1 in the recent work by 
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McOuinn et alj (|2007|) . Infrared Array Camera (IRAC) magnitudes, obtained in five different 
epochs spanning 19 months from January 2004 to January 2006, in the 3.6 and 4.5 /im channels 
can be found in this last reference; UIT 005 was not detected at 8.0 /im. The values shown in Tab|T] 
are the mean values of those five individual measurements. The uncertainties in the individual 
me asurements are compar able to the error in the mean values, which is consistent with the fact 



that 



McQuinn et al 



(|2007r) did not consider UIT 005 among the variable stars in their sample. 



The photometric measurements, along with the corresponding dates, are summarized in 
Tab.[IJ We show in Fig. [Da false color image of a 1.5'xl.5' field approximately centered 
on UIT 005, iising B (blue), V (green) and Ha (red) images from the Local Group survey 



(IMassev et al 



2006). Fro m the low level variations in V and B-V, and the non-variability assigned 



by lMcOuinn et al.l (120071) . it seems reasonable to conclude that the star is not showing photometric 
variations beyond the stochastic level. The changes in the U-band would be most likely related 
to issues with the calibration of the different ground-based measurements. As it will be shown 
later (see all these photometric measurements, excluding the U-band, are well reproduced 
by the spectral energy distribution (SED) predicted by the model corresponding to our derived 
fundamental parameters. 



2.2. Spectroscopy 

The quantitative analysis presented in the following sections is mainly based on two 
observational runs with the 4.2m William Herschel Telescope (WHT) equipped with the 
Intermediate dispersion System (ISIS). The blue wavelength range (3970-4850 A) was acquired 
in September 7th 2000, using a 1200 lines mm~^ grating (dispersion, 0.23 Apix^^) centered 
at 4400 A. The 1.2 arcsec slit was imaged with an EEV12 detector (scale, 0.19 arcsecpix ^), 
providing a spectral resolution of ~0.9 A full width half maximum (FWHM). Two 3600 s 
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exposures were secured under a seeing of 2 arcsec, for an airmass below 1.3. 

The red spectrum was observed in October 15th 2003, using ISIS red arm equipped with 
a 600 lines mm~^ grating (dispersion, 0.44 A pix-i). A 1.2 arcsec slit was used, imaged with 
the MARC0NI2 chip (scale 0.20 arcsec pix^^), providing the system a spectral resolution of 
1.6 AFWHM. Two exposures, one of 2200 s and another of 2000 s, were acquired, for a char- 
acteristic airmass below 1 .4. The red spectrum expands approximately 1 200 A, centered at 6558 A. 



A third data set was collected from the Keck Observatory Archive. The star was observ ed 



in Oc tober 10th, 2001 with the High Resolution Echelle Spectrometer (HIRES, 



Vogt et al. 



19941) at Keck 1. This set contains 4 exposures of 3600 s each, covering the range ^4265-6630 
A at R~65000, with gaps in the coverage between some of the 30 orders. Unfortunately, the 
instrumental set-up used during these observations was not optimal for B-type stars, and a number 
of key diagnostic lines, in particular the Si II and Si IV lines that are required to determine the 
effective temperature, were not included. Whilst an independent analysis of this spectrum is then 
not possible, we can still compare it with the ISIS spectra, to search for spectral variations. 



The ISIS data were processed using standard IRAFlfl routines for long slit spectroscopy (bias 
correction, flat-fielding, wavelength calibration and optimal extraction). The archival HIRES 
data were fully processed with Tom Barlow's MAKEE package. The final extracted spectra were 
continuum normalized using low-orders polynomial fits to well defined continuum points. 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by 
the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



New, service mode, observations with ISISAVHT were carried out on November 6th, 2009, 
under clear sky conditions, for a seeing of ~1.2 arcsec and a characteristic airmass of 1.7. In 
this occasion, both arms were fed simultaneously. The configuration of the blue arm was quite 
similar to the one used in the 2000 run, but with the spectrum centered at 4260 A, hence covering 
the range 3868-4618 A. For the red arm, we selected the R1200R grating, centered at 6800 A, 
with the spectrum covering the range 6497-7190 A at a spectral resolution of ~1 A FWHM for 
a projected slit width of 1.2 arcsec. Two exposures of 1800 s each were collected. This set was 
processed in the same way as explained above. The S/N achieved (S/N~30-40) is well below the 
one in the 2000/2003 cases (S/N~80) ; nevertheless, this new spectrum shows that the star has not 
changed. 

We have carefully compared the spectra from the various epochs for radial velocity variations, 
finding no significant evidence for binarity. However due to difficulties in assigning zero-point 
off-sets for the various observational configurations we are unable to detect radial velocity 
variations smaller than ±10 kms^^ . 



3. Spectral Analysis 

3.1. Qualitative description of the spectrum 

The first noteworthy spectral characteristic is the narrow H Balmer profiles, pointing to a 
high luminosity. This is supported by the 2003 Ha profile, fully in emission. In terms of metal 
features, the blue spectrum is characterized by the presence of strong N II lines, and the weakness 
of the O II lines. Other prominent features present in the spectrum are the strong Si III triplet 
lines at 4553-67-74 A and the Mg II doublet at 4481 A. Three other strong N II lines appear in 
the observed ISIS red domain, along with Ha and He I 6678 A line, plus the C II doublet at 



- 11 - 



6578-82 A. 

The two Ha profiles, shown in Fig. |2l are slightly different. In 2001, the profile (top panel in 
Fig. [21) appears as a P Cygni line, whilst the 2003 observation (mid panel) shows a strengthening 
of the profile in pure emission. To better illustrate this change, we have included in the bottom 
panel of Fig. [2] a comparison of both observations, with the HIRES 2001 profile degraded to the 
spectral resolution of the ISIS 2003 observation. By contrast, the few metal lines in the the red 
range do not show any change, within the noise level. In the same sense, there are no noticeable 
differences in the blue range between the ISIS 2000 and HIRES 2001 spectra. 



By comparing the equivalent wid th ratios of Mg II 44 8 1 A to the Si III triplet li nes with the 



Fitzpatrick 



sampl e of Galactic supergiant stars of Lennon et al.l (|I993l) and the LMC sample of [ 
(1199 lb . we tentatively locate the star in the spectral type range of B2-2.5. Although the metallicity 
of UIT005 could be different from th at of these templates , the abundance ratio of these two 
elements remains fairly constant (see 



Urbaneja et al 



2005|) . Correspondingly, we assign a 



luminosity class of Ia+ due to the narrow H Balmer lines. Thi s is a minor revision in luminosity 



Massev et al. 



(119961) on the basis 



class with respect to the classification previously provided by 
of low resolution spectroscopy (see their Fig 5b, displaying a spectrum covering the wavelength 
range ~3900-4900 A, obtained either in 1993 or 1995). 



3.2. Methodology: model atmosphere/line formation code and main assumptions 



The model atmosphere/line formation code FASTWIND (Fast Analysis of STellar atmospheres 
with WINDs) was used to create a small set of mod els for the analysis. The code was first 



introduced by 



Santolava-Rev. Puis & Herrero 



(|l997|) . Meanwhile, a number of improvements 



have been incorporated. A complete description of the current status of the code as well as 
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comparisons with alternative codes have been presented by IPuIs et alj (120051) . Briefly, as used 
in this work, FASTWIND solves the radiative transfer problem in the co-moving frame of the 
expanding atmospheres of early-type stars in a spherically symmetric geometry, subject to the 
constraints of statistical equilibrium and energy conservation. Steady state and homogeneous 
chemical composition are also assumed. The density stratification is set by the momentum 
equation (the advection term is not considered) along with the mass-loss rate and the wind 
velocity field (a standard /3-law) via the equation of continuity. A smooth transition between the 
wind regime and the pseudo-static photosphere, which is described by a depth-dependent pressure 
scale-height, is ensured. 

Each model is defined by a set of parameters: effective temperature, surface gravity and 
stellar radius (all these quantities are defined at Tross = 2/3), the exponent of the wind velocity 
law (/3), the microturbulence velocity, the mass-loss rate, the wind terminal velocity and a set 
of chemical abundances. In order to reduce the num ber of parameters to b e explored with our 



Urbanej a et al 



(120051) . with a modification 



grid of models, we follow the procedure described in 
in that the microturbulence velocity is fixed to 10 kms ^ for the calculation of the atmospheric 
structure. For specific information concerning the model atoms used in our calculations, the 
reader is referred to the last reference. 



3.3. Stellar parameters and surface chemical composition 

Tab. |2] summarizes the derived properties of UIT005, with Fig. [3] and Fig. |4] displaying 
respectively a comparison of the ISIS and HIRES data with our final model. The effective 
temperature is well constrained by using the Si II / Si III ionization balance, and the surface 
gravity is obtained from the fits of H5 and Hy. The shape of the emission Ha profile allows us 
to constrain the wind- velocity law (/3), although the profile is not completely well reproduced, 
in particular the blue ward absorption. This points to limitations in some of the assumptions in 



our modeling technique, in particular related to the density distribution. However, the derived 
mass-loss rate is not seriously affected by these possible deficiencies due to the fact that the 
emission part of the profile reacts strongly to changes in this quantity. Also, both Hy and H5 
profiles appear well reproduced, which reinforces our confidence on the global procedure. 

In order to constrain the He abundance, we consider the He I lines at 4026, 4713 A (both 
belonging to the triplet spin system) and the one at 4387 A (singlet spin system), while 6678 and 
4471 A (singlet and triplet systems, respectively), although observed, are not directly considered. 
The lines at 4009 (singlet), 4129 (triplet), 4145 and 4437 A (both singlet lines) are also not 
considered because presently our models do not include detailed data for the Stark broadening 
profiles for these lines. 

Concerning the N abundance, lines from single ionized N are used to determine its surface 
abundance. These are 3995, 4447, 4601-07-14-21-30 A along with the lines present in the red 
wavelength coverage, namely 6379, 6482 and 6610 A. Other N II lines can be identified at 
4035-41-43 A, but they involve high energy terms, presently above the energy cut-off considered 
in our N II model. 

The derived carbon abundance is based on the three available lines: C II 4267, 6578-82 A. It 
must be noted that we regard our current C model atom as not fully reliable, and thus the quoted 
C abundance has to be considered carefully. 

For the determination of the O abundance we consider the single ionized oxygen lines located 
at 4072-76, 4317-19-66, 4414, 4591-96, and 4661 A, neither of them known to present any blend 
with any other species in the range of effective temperatures characteristic of our star. Other O II 
lines, not directly used, are present in the spectrum, such as the lines around ~ 4650 A, the weak 
lines at 4673-6 A and, particularly, the group of strong lines just redwards of the red wing of Hy 
(located at 4345, 4347, 4349 and 4351 A). In general this group of lines forms on the red wing of 
Hy, whilst due to the narrow H Balmer profiles, these lines appear out of the wing in this case. 
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Finally, the Mg abundance is based on the only one feature present in the spectrum, the Mg II 
line at 4481 A. 



3.4. Error analysis 



Uncertainties affecting the determination of the effective temperature and the (effective) 
surface gravity basically depend on the quality of the observed spectrum (i.e. the signal-to-noise 
ratio). In our case, Teff can be constrained within a range of ± 1000 K : any larger variation in 
this quantity would shift the Si ionization equilibrium in such a way that either Si IV or Si II lines 
would became too strong relative to the observed ones. 

The effective surface gravity is well determined by simultaneously fitting the H Balmer lines 
(H5 and Hy in this case), from which we found an uncertainty of ±0.05 dex. To that, we have 
to add the contributions coming from the stellar wind parameters. In doing so, the total budget 
amounts to ±0.10 dex, taking into account ver y conservative errors for /3 and Q' (the optical 



depth invariant, see f.e. 



Kudritzki & Puis 



2000[ is defined as Q' = M/J/^ {R^vJ)^^'^ , where 



/ci represents the clumping factoci set to unity in our models). We would like to stress out that 
here we are disregarding any possible contribution from the models themselves, i.e. we consider 
that our models properly account for the photospheric line pressure. This assumption is related 
to the com pleteness of the line list used in the blanketing/blocking calculations (the reader is 



referred to 



Puis et al, 



2005 



for a complete discussion concerning this point). A second issue is 
related to the lack of information about the Fe abundance. Here we are considering that the Fe 
abundance (actually the abundance of any element other than He, C, N, O, Mg and Si) scales with 



^Presently, clumping is treated under the micro-clumping formalism: the wind consists of over- 
dense clumps embedded in a void inter-clump medium; /d measures the overdensity in the clumps. 
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the a-elements according to the solar ratio and the global metallicity of the model. This could 
have an impact on the effective gravity (again, through the contribution to the photospheric line 
pressure). As a sanity check, we have tested that varying the Fe-group element abundances by 
±0.2 dex, no noticeable effect is observed on the synthetic H profiles used for the determination 
of the surface gravity. 

Finding the stellar radiu s requires the knowled ge of the absolute magnitude, for which the 



distance to the star (840 kpc. 



Freedman et al, 



200 ih and photometric data in several bands are 



required. Along with My, the radius depends also on the SED provided by the models, and in 
particular on the flux in the Johnson V-band. The error in the absolute magnitude is dominated in 
our particular case by the uncertainty of the distance modulus to the galaxy, A/i = ±0.10 mag 
(see previous reference), since in comparison the adopted error in the apparent magnitude is 
small. Amy = ±0.05 mag. This AMy = 0.10 mag translates into a contribution of AR(/i) = 4 
Rq due to the uncertainty in the distance. We can also estimate the contribution due to the 
differences in the theoretical SEDs (basically the uncertainties on Teff) by computing the required 
radii to reproduce the absolute magnitude for our limiting models with Teff± ATeff, resulting in 
AR(Teff) = 4R0 . Considering both contributions, we estimate an uncertainty in the radius of 
AR = 6R0. 



In principle, it would seem feasible to estimate the wind terminal velocity from the Ha P 
Cygni profile present in the HIRES spectrum. In doing so, the velocity at which the continuum 
level is recovered corresponds to ~125 kms^^ . However, the P Cygni profile is not saturated, 
therefore it is unclear whether this is the actual terminal velocity; it could happen that the opacity 
in the line beyond this point is negligible and it is not contributing to the formation of the line. 
Models computed for a broad range of terminal velocities (keeping Mv^^ constant) produced 
indistinguishable profiles, hence suggesting that the line should not be used to estimate this 
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quantity. Without precise information of the wind terminal velocity, the goodness of our derived 
mass-loss rate M, is unsure. It is possible to assign formal errors to M by model comparisons, 
but it is important to understand that the derived mass-loss rate is anchored to the adopted wind 
terminal velocity, in the sense that two models with different M and Voo but keeping Mv^^ constant 
will show very similar Ha profiles. 



2002 



2000 



Urbane] a et al 



Vink et al 



2002 



20001 : 



Evans et al, 



Kudritzki 



20021 : 



R^ecent UV studies of B-type supergiants in different metallicity environments (iBresolin et al. 



2004) aswel 



Krticka 



as theoretical investigations (iPuls et al 



200a), indicate a weak metallicity dependence 



of the terminal velocit y. On the other hand, observational studies of Galactic stars (see 



Kudritzki & Puis 



2000l for a review) show a relatively large dispersion in wind terminal velocities 



within a given spectral type. For the sake of the analysis, we have adopted a terminal velocity 



given 



by the surface escape velocity and the Vcx,-Vesc relationship provided by 



Kudritzki & Puis 



(|2000|) . The mass-loss rates quoted in Tab. |2] correspond to the values obtained from the two 
different Ha profiles available, and, to avoid any over interpretation, we do not include a precise 
estimation of their uncertainties. 



Finally, at the spectral resolution used in this work, the signal-to-noise ratio is the main 
contributor to the uncertainties in the chemical abundances. Other factors could increase the 
line-to-line scatter for each element, like the continuum normalization and the accuracy of the 
derived fundamental stellar parameters. The uncertainties presented in Tab. |2l result from the 
combination (quadratic sum) of two different contributions. The first one corresponds to the 
statistical error associated with the elemental abundances derived from our final model, based on 
a line-by-line analysis (see Tab. [3]). For the second one, we consider the uncertainties related to 
the determination of the fundamental parameters, taking into account the uncertainties affecting 
the effective temperature, the surface gravity (which included the effect of wind parameters), the 
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microturbulence and the uncertainty in the helium abundance. These individual contributions are 
compiled in Tab. H where the 1 o uncertainties quoted in the last column of the table are obtained 
as the quadratic sum of the corresponding contributions. 



4. Discussion 

Besides the faint source already discussed in Sect. 12.11 there is no indication in the different 
spectra that would suggest that UIT 005 is in fact a multiple system. If indeed the star has 
companions, these are not contributing at any significant level. There is also no reason to assume 
that our object and the nearby faint source are physically related in any way. In the following, we 
consider that UIT 005 has evolved as a single star, since this is the simplest assumption that can be 
made from the information in hand. Moreover, we adopt the characteristic metallicity of UIT 005 
as given by the Si and Mg abundances, ~0.7 Z©, i.e. a metallicity intermediate between that of 
the Milky Way and that of the LMC. This value, based on the abundances of two a-elements, 
is consistent with what would be inferred from the C+N+0 by mass relative to the solar value, 
~0.8 Zq. We note here that this metallicity is higher than the corresponding value for the ionized 
gas in NGC 588 when using the gas phase O abundance (see below) as a proxy for metallicity. 



OAZgy . The population synthesis based study of the ionizing cluster of NGC 588 by 



by 



Jamet et al, 



(120041) indicates an eve n lower metallicity, O.TLq. However, both these results are contradicted 



Bianchi et al. 



(|2004r) . who find near solar metallicity for UIT 008, one of the two known 



Wolf-Rayet stars in NGC 588. 



Fig. [5] locates UIT 005 in a log Teff-log g di agram, together with the sample of M33 blue 
supergiant st ars analyzed by Urbaneia et al \2OO5 , hereinafter UHK05), and evolutionary models 
for Galactic ([Meynet & Maedeil 120031) and LMC ([Meynet & Maedeil 120051) metallicities. The 



position of UIT005 is consistent with models having initial masses between 40 and 60 for 
both metallicities. As the derived metallicity of UIT005 lies between that of the Galaxy and the 
LMC, we conclude that its initial mass also was between these two values. Note that, unlike in 
the classic Hertzsprung-Russell (H-R) diagram (luminosity versus effective temperature), both 
Teff and logg are directly derived from the spectral analysis, without any assumption regarding 
the distance to the object. This figure also illustrates that UIT 005 is farther away from the Main 
Sequence phase than the M33 B-type supergiant stars in the sample studied by UHK05. This 
point is stressed in Fig.|6l our object appears clearly isolated in this N/O-Teff diagram, showing a 
much higher N/0 ratio (by mass) than any of the other M33 objects. The star presents a N/0~8, 
and a N/C~32 (this last value is highly uncertain because of the issues with our C model atom). In 



compa rison, the object with the higher N/0 ratio in the Galactic sample studied by 



Crowther et al. 



(120061 empty circles in Fig. |6]) is the well known blue hypergiant HD 152236, for which these 
authors find N/0 = 3. As we will discuss later, the apparent agreement with the predictions by the 
evolutionary models is misleading. 

Even though the abundances are far from the expected values at (CNO) equilibrium, 20.5 
and 64 respectively, the very high observed N/0 ratio of UIT 005 is quite remarkable and, at a 
first glance, could be an indication of an advanced evolutionary status (enhanced N abundance). 
Altematively, this could be a consequence of an initial low O abundance, be cause of the lo cation 

and its strong radial abundance gradient (UHK05; 



of the star in the galaxjc 
could result from a combination of both. 



U et al. 



2009 ). Or it 



The proposition of an initial low O content can be easily tested by using the information 
provided by the other two a-elements for which we derived an abundance, Mg and Si. Fig. |7] 
displays the O and Mg abundances of UIT 005 and the sample studied by UHK05. In the case 



'^Normalized distance to the galactic center, R / R25 = 0.80 
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of O (Fig.|71 top panel), it seems that UIT005 does not follow the gradient defined by the other 
B-type supergiants, whilst this is not the case for Mg (Fig.|7l bottom panel). The same happens 
if we consider Si instead of Mg (not shown). Using the stellar O radial gradient obtained by 
UHK05, the O abundance that would correspond to the galactocentric distance of UIT005 is 
log(0/H)+12 = 8.45 dex. Even if one argues that there is no particular reason for the gradient to 
behave smoothly, our expectation for the O abundance is supported by recent results for the giant 
H II region NGC 588, separated by ~30" (a de-projected distance of ~130pc for the distance 
to the galaxy adopted in this work) of our star, with an O abundance in the range 8.45-8.32 dex 



dCrockett et al 



2006 : 



Magrini et al. 



20071 



Rosolowskv & Simon 



2008). It would seem then that 



the present low O content of UIT 005 is not related to its location in M 33. 

We can further investigate this issue by removing the spatial dependence. Fig. [8] presents 
the relation between the derived O abundance and the mean a-element abundance defined by 
Mg and Si. Considering the sample of UHK05, it seems clear that there is a relation between the 
abundances of these three elements, that is not followed by the abundances of UIT 005. Therefore, 
the relative stellar abundance of O to the two other a-elements, Mg and Si, in combination with 
the previous paragraph, strongly support the idea of UIT 005 having suffered significant evolution 
and presently exposing material processed by the CNO cycle. 



Further constraints can be obtained from the evolutionary lifetimes of phases consistent 
with the effective temperature and gravity derived for UIT 005. Not all the possible models have 
the same likelihood, since in some cases the compatible stages are predicted to be extremely 
short. In fact, a close examination of the lifetimes indicate that for the rotating models these are 
always only a few hundred years, making it extremely improbable that UIT 005 started its life 
with a high rotational velocity. On the other hand, non-rotating models at 50 and 60 Mq spend 
respectively 6.7 and 6.2 x 10* years in this phase at Galactic metallicity, while the 50 M© at LMC 
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metallicity spends 3.8x10'^ years, giving us a chance to observe such an object. Therefore, from 
the available set of evolutionary models, we are left with the non-rotating models at 50 Mq (for 
both metallicities) and 60 Mq (only for the solar metallicity, but note that the non-rotating 60 M© 
model was not available for LMC metallicity). For all three models the evolutionary predicted 
masses are between 30 and 44 M©, in agreement with the spectroscopic mass, within the error 
bars. At the same time, however, the observed surface enrichment is at odds with the ones 
provided by these evolutionary models; whilst N is within the uncertainties of the observed value, 
the models predict that the O surface abundance should be almost an order of magnitude higher 
than observed. Since, owing to rotation, one would expect a higher degree of evolution in N than 
in O, the observed values would not be explained by models with higher initial rotational velocity 
alone. 

Because of the somewhat peculiar derived abundances, one might ask to what extent our 
results are conditioned by the necessary assumption incorporated in the model atmosphere/line 
formation code used for the analysis. During th e course of our work, we calculated a number 



of models with the altemative code CMFGEN (Hillier & Miller 



1998 ). These models were 



based on our fastwind solution for UIT005; i.e. we did not perform an independent analysis 
using CMFGEN models in the same way as described in previous sections. We limited the 
calculations to a small region of the parameter space around the values obtained with FASTWIND 
instead. A discussion of the differences between both codes is clearly beyond the scope of 
our paper; this would require a larger sample of objects, to investigate possible systematic 
differences. Nevertheless, the "best" solution found with CMFGEN (Teff ~18300 K, logg ~2.15 
dex) overlaps, within the errors of the analysis, with our parameters. Moreover, the chemical 
abundances derived with both codes are consistent. In particular, the O abundance obtained 
with CMFGEN, 7.85 dex, is in excellent agreement with our fastwind derived value of 7.90 dex . 
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Although somehow uncertain due to the lack of precise information about the wind terminal 
velocity, the mass-loss rates derived from the available Ha profiles (M ^ 1.0 — 1.4 x 10~ ^ 



Vink et al, 



(12000 



Ma,y r^^ ) are low when compared with the theoretical value predicted by 
200 lb . M = 1.8 X 10^^ MQyr^^ . We have not found firm indications of a structured wind in 
our admittedly limited set of observations. But even if the wind is clumped to some extent, this 
will not bring the observed and the theoretical value into a better agreement. Rather, the situation 
would become even worse, since a clumped wind would act in the opposite way, implying a 
smaller true mass-loss (our derived mass-loss rates would be upper limits in that case). This is a 
significant poi nt of disagreement with the evolutionary models, since these apply the theoretical 



predictions by 



Vink et al. 



(2000 



200 lb to prescribe the mass-loss rates during the blue supergiant 



phases. This discrepancy between theoretical predictions and Ha and radio continuum derived 



mass-loss rat es for mid-B types has been previously reported for SMC (ITrundle & Lennon 



and Galactic dCrowther et al 



refer the reader to 



2006 



Markova & Puis 



1 r 

Markova & Pulsli2008l: iBenaglia et al 



2005 ) 



2008b supergiants. We 



(|2008|) . where this dilemma has been discussed at some length. 
In fact, in the case of UIT 005, the high theoretical mass-loss rate is not consistent either with the 
observed spectral energy distribution (see below). 

The available Ha profiles indicate a variation of ~25% in the mass-loss rate from the 2001 
to the 2003 observations. At the same time, the few He and metal lines available in the same 
wavelength range suggest that this change does not affect the structure of the star significantly. 
Changes in Ha profiles at this level have been previo usly reported for Galactic early B-type 



supergiants ([Rivinius et al. 



19971: 



Crowther et al. 



2006 ). 



UIT 005 resides in a region of the H-R diagram (Fig.|9l) populated by extreme blue objects 
such as LB Vs. In fact, if the dista nce to M33 is larger than assumed here, as recently derived 
from several stellar indicators (see iBonanos et al.ll2006t lU et al.ll2009l and references therein). 
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then UIT005 could be as luminous as log(L/L0) ~ 6.1dex, well in the regime of the most 
luminous LB Vs. Whether our star is a regular blue object or a dormant LBV is unclear. The 
observed spectral energy distribution (SED) from 0.44 to 4.5 /im (Fig.fTOl) is very well matched 
by the synthetic SED of our final model (Fig.fTOl solid line), even though those photometric 
measurements were acquired in different epochs. The synthetic SED has been reddened 
assuming a total-to-selectiye extin ction ratio Ry = 3.1 and by adopting the extiction curve by 



Cardelli. Clavton & Mathis 



(|1989l) . We note that the selection of the extinction curve has little or 
no effect on the considered photometric bands, and the excellent agreement between observed 
and theoretical SEDs would not be affected. We have included a second model (dash-dotted line 
in Fig.fTOl) computed for an enhanced mass-loss rate, set to the theoretical prediction discussed 
above. As previously indicated, the observed SED for wavelengths longward the Lband is clearly 
not consistent with this high mass-loss rate. 

With respect to our final model, the good agreement in the infrared bands, and in particular 
the no detections at 8.0 and 24 /im (this last from inspection of archival Spitzer data), would 
indicate the lack of a significant amount of hot dust around the star, suggesting that UIT 005 is not 
surrounded by a nebula. Whilst the presence of such nebula is a common characteristic for many 
LBVs and it would support the LBV nature of UIT 005 , there are also cases of confirmed LBVs 



(for example Rl 10, Sk -69° 142a, and R85 in the LMC, 
absent. 



Bonanos et al 



20091) where this nebula is 



The nature of this intriguing object is further complicated by its apparent isolation, as 
illustrated by Fig. [T] where we see that it lies well beyond the periphery of its nearest cluster, 
NGC 588. Is it possible that this cluster was the original birth-place of UIT 005? In Fig. [TI] 
we show the color-magnitude diagrams (CMD) for all stars in the fie ld (aproximately cov ering 



2.7x2.7 arcmin^). It is known that NGC588 contains two WN stars (|Drissen et al 



2008) so it 
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is potentially both young enough and massive enough to have produced a supergiant such as 
UIT 005. As Fig. [m demonstrates, UIT005 fits nicely into the CMD of NGC 588 as one of the 
most evolved of its potential members. Furthermore, ass uming a main-sequ enc e lifetime of aroun d 



5 My r (consistent with the age of 4.2 Myr estimated by 



Jamet et al. 



2004 and 



Ubeda & Drissen 



2009L for a somewhat lower metallicity), it would require a runaway velocity of around 40-50 
kms^^ to reach its current dist ance from the cluster, well within the range of velocities of known 



Galactic runaway 0-type stars ([McSwain et al, 



20071). While not conclusive, for example there 



are no radial velocity measurements for stars in NGC 588, the evidence above provides interesting 
circumstantial support for the idea that UIT 005 is a runaway star from NGC 588, with the ejection 
mechanism being either dynamical interaction in the cluster, or a SN kick (from a yet more 
massive companion). This second scenario would require a much higher runaway velocity, since 
the SN progenitor would need some time to evolve. Whilst this seems to favor the dynamical 
interaction as the ejection method, very little is know about the cluster itself, hence this scenario 
cannot be investigated further at this point. 



5. Conclusions 

Our detailed quantitative analysis of UIT 005 shows that this massive and luminous blue 
supergiant (B2-2.5Ia+, logL/L© ~ 5.9dex, and M~5OM0) has a surface abundance pattern that 
presents peculiarities, in particular when considering its location in the disk of the M 33. From the 
point of view of the a-elements for which we have information, Mg and Si abundances seem to 
follow the radial behavior (as traced by other B-type supergiants), while its O abundance appears 
to be depleted by a factor ~3.5 with respect to the expected local value. The a/0 ratios are 
hardly understandable as a product of previous stellar generations at this particular location in the 
galaxy, moreover if we incorporate the derived high N/0 (~8, by mass) and the He abundance. 
All together, and under the assumption that the star is evolving as a single object, the most likely 
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explanation is that the observed chemical pattern is a consequence of the chemical evolution 
of the star, that is showing an advance stage of CNO processing in its surface. This picture is 
supported by the predictions of recent evolutionary models, although we must point out that some 
discrepancies are still present, such as the degree of chemical processing predicted by the models. 

The derived parameters place the star in a region of the H-R where LBVs are usually found. 
The observations do not support its nature as an LBV, nor they rule out the possibility of an LBV 
in a dormant state. UIT 005 could very well represent a transitory stage in the evolution of the 
most massive stars or even the precursor of an LBV. Subsequent monitoring of this object would 
help us to improve our knowledge of massive stars, bridging the gap between regular and more 
exotic blue supergiants. 
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Table 1 . Photometric data 



Magnitude 



Bandpass/Color 


(mag) 


Date 


Comment 


V 


16.78 


1993-11-22 


Massev et al. (1996), around 


B-V 


0.03 


1993-11-22 


Massev et al. (1996). ground 


U-B 


-0.91 


1993-11-22 


Massev et al. (1996). around 


F170W 


15.566±0.013 


1994-10-28 


this work, WFPC2/HST, no transformation 


F336W 


15.415 ±0.004 


1994-10-28 


this work, WFPC2/HST, U= 15.970 ±0.004 


F439W 


16.813±0.004 


1994-10-28 


this work, WFPC2/HST, B=16.833±0.004 


F469N 


16.845±0.014 


1994-10-28 


this work, WFPC2/HST, no transformation 


F547M 


16.836±0.003 


1994-10-28 


this work, WFPC2/HST, V=16.831 ±0.003 


V 


16.727 ±0.003 


2001-09-18 


Massev et al. (2006), ground 


U-B 


-1.053 ±0.003 


2001-09-18 


Massev et al. (2006), around 


B-V 


0.010±0.003 


2001-09-18 


Massev et al. (2006), ground 


V-R 


0.019±0.003 


2001-09-18 


Massev et al. (2006), ground 


R-I 


0.058 ±0.003 


2001-09-18 


Massev et al. (2006), ground 


J 


16.79±0.135 


1997-12-05 


Cutri et al. (2003), 2MASS 


J 


16.77±0.01 


2005-09-29-2005-12-16 


Cioni et al. (2009). WFCAM/UKIRT 


H 


16.74±0.01 


2005-09-29-2005-12-16 


Cioni et al. (2009). WFCAM/UKIRT 


Ks 


16.77±0.03 


2005-09-29-2005-12-16 


Cioni et al. (2009). WFCAM/UKIRT 


IRAC 3.6;im 


16.637 ±0.078 


2004-01-09 to 2005-08-25 


McOuinn et al. (20071 mean values 


IRAC 4.5 


16.659±0.145 


2004-01-09 to 2005-08-25 


McOuinn et al. (2007). mean values 
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Table 2. Stellar parameters and chemical abundances 



Properties UIT005 



Alt. ID 


J013242.92+303847.0 


SpT 


B2-2.5 Ia+ 


Teff (kK) 


19.0±1.0 


log g (cgs) 


2.25±0.10 


Voo (km ' ) 


450 




1.42-1.00 


R(R0) 


86.5 ±6.0 


P 


2.20 ±0.20 


n(He)/n(H) 


0.15 


v»„-A (kms-') 


10±3 


V sin; (kms^' ) 


< 50 


A(C) 


7.40 ±0.20 


A(N) 


8.85±0.15 


A(0) 


7.90 ±0.22 


A(Mg) 


7.50±0.16 


A(Si) 


7.35±0.15 


Z(Z,;,)=' 


0.75 


m,, (mag) 


16.79±0.05 


B - V (mag) 


0.01 ±0.01 


E(B-V)(mag) 


0.18±O.O15 


Mv (mag) 


-8.33±0.11 


EC (mag) 


-1.79±0.12 


Log (L/L,,) (cgs) 


5.95±0.16 


M'^P'='-' (M0) 


48.6±18.0 



"The chai'acteristic metallicity Z is defined as 
the mean of the individual differences of Mg and 
Si abudances with respect to their solar values 

Note. — The mass-loss rates correspond to the 
values derived from each Ha profile. 
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Table 3. Equivalent widths, line-to-line abundances and statistical uncertainty of the mean 

values 





EW 






Fcsturc 


fmA'l 


SNR 


(dex) 


Cll 4267 


104 


60 


7.37 


Cll 6578 


154 


70 


7.44 


C7(C) 






0.07 


Nil 3995 


365 


34 


9.06 


Nil 4447 


117 


57 


8.25 


Nil 4601 


177 


65 


8.75 


Nil 4607 


173 


62 


8.85 


Nil 4614 


144 


60 


8.82 


Nil 4621 


189 


86 


8.91 


Nil 4630 


358 


67 


9.00 


Nil 6380 


123 


84 


8.75 


Nil 6482 


316 


63 


8.87 


Nil 6610 


140 


99 


8.51 


CT(N) 






0.08 


On 4072 


48 


50 


7.84 


Oil 4076 


87 


50 


7.54 


On 4317 


32 


60 


7.63 


On 4319 


93 


60 


8.40 


On 4366 


42 


70 


7.76 


On 4414 


81 


65 


7.97 


On 4591 


22 


55 


7.38 


On 4596 


22 


50 


7.53 


Oil 4661 


100 


58 


7.95 


a(0) 






0.11 


Mgll4481 


176 


74 


7.50 


CT(Mg) 








Sim 4552 


272 


68 


7.32 


Sim 4567 


217 


75 


7.31 


Sim 4574 


155 


80 


7.38 


CT(Si) 






0.04 
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Note. — Abundances derived for each 
line; SNR as measured in tlie continuum 
close to the hne. 
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Table 4. Metal abundance uncertainties: contribution from different sources 





+(Ar,Alogg) 


-(Ar.AIogg) 


+ Avturb 


- Avturb 


-AY 


-A log Q' 


a 


Species 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


(dex) 


C 


0.09 


-0.16 


-0.02 


0.02 


-0.04 


-0.10 


0.19 


N 


-0.07 


0.05 


-0.10 


0.09 


0.00 


-0.03 


0.13 





0.05 


-0.17 


-0.04 


0.04 


-0.03 


-0.06 


0.19 


Mg 


-0.10 


0.11 


-0.06 


0.06 


0.02 


-0.09 


0.16 


Si 


-0.03 


-0.14 


-0.03 


0.04 


-0.02 


0.02 


0.15 



Note. — The last column presents the 1 (J uncertainties obtained as the quadratic sum of the different 
contributions. These Iff uncertainties are combined with the line-to-line abundance scatter (see previous 
table) to produce the finally adopted metal abundance uncertainties quoted in Tab.|2] 



Fig. 1. — (Left) False color composite image of a 1.5 x 1.5 arcmin^ field around UIT005. North 
is up and East is to the left; UIT005 is the blue bright source slightly to the West of the image's 



center. The nearby giant H II region NGC 588 is inclu ded. This view resul ts from the combination 



of three images from the Local Group Galaxy Survey (iMassey et al, 



20061) : red-Ha, green- V, and 



blue-B. (Right) A wider field is shown with the WFPC2 footprint overlaid on a V band image; 
UIT 005 is identified by a red cross. The orientation of both images is the same. 
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Fig. 2.— UIT005 Ha profiles. From top to bottom: HIRES (2001, R~65000), ISIS (2003, 
R~5000), and a comparison of the two observations with the HIRES profile (thick) degraded to 
the ISIS resolution (thin profile). 
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Fig. 3. — Observed ISISAVHT spectrum (thin line) of UIT 105 and the synthetic model (thick line) 
resulting from our analysis. An identification of some prominent features is provided. 
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Fig. 4. — Observed HIRES/Keck 1 spectrum (thin line, R~65000) of UIT005 compared with our 
final model. Only selected windows are shown, to illustrate the excellent agreement between the 
model based on the lower resolution ISISAVHT spectra and the high resolution HIRES spectrum. 



-39- 




Fig. 5. — Derived Teff-log g c ompared with the predicti ons of stellar evolutionary models ac 



counting for stellar rotation from 



Meynet & Maederl (|2003|) for solar metallicity 



masse s 9, 12, 15, 20, 25, 40 and 60 M© are displayed as solid lines) and from iMeynet & Maeder 



(tracks for ZAMS 



(l2005h for z=0.008 (30 and 40 Mq tracks are sho wn as dashed-dotted li nes, and are labeled at the 
upper right comer). The B-type supergiant stars of lurbanej a et aJ ( 2005 ) are shown as circles, with 
UIT 005 represented by the square. 
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Fig. 6. — Surface N/O ratio (by mass) against the effective temperature. Theoretical predictions 
and filled symbols as in previous figures. For reference, the sample of Galactic B-type supergiants 
studied by ICrowther et al.l (|2006h is also shown (small symbols without error bars). 




Fig. 7. — UIT005 oxygen and magnesium abundances in the context of M33 radial metallic- 
ity gradient. B supergiant stars from UHK05 are shown as circles, with the square representing 
UIT005. The x-axis corresponds to the normalized de-projected distance to the galaxy center, for 
the distance to M33 and its structural parameters defined in fJT] The 1 o confidence intervals for 
the gradients are shown by the thin lines. 
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Fig. 8. — Oxygen vs magnesium plus silicon abundances, relative to the Sun. The square represents 
the abundances of UIT 005, and the circles correspond to the sample analyzed by UHK05. A linear 
fit to this sample is shown by the dashed line; the la confidence intervals are shown by the thin 



lines. Solar abundances from 
O, Mg and Si, respectively. 



Asplund et al, 



(l2005h . namely: 8.67 dex, 7.56 dex and 7.54 dex for 
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Fig. 9. — Hertzsprung-Russel diagram. The square represents UIT005, whilst the big dots locate 
the sample of early B-type supergiants studied by UHK05. The approximated location of the S Dor 
instability strip and the maximum light strip are also indicated. Known LBVs and LBV candidates 
are included in this plot (small symbols). When the information is available, the two different 
states of a given object are connected by dots. Evolutionary tracks for Z = 0.008 from the Geneva 
group are displayed, for different initial m asses and an initial equa t orial r otation of 300 kms^ . 



This figure has been partially adapted from 



Smith. Vink & de Koteil (l2004h . 
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Fig. 10. — UITOOS's spectral energy distribution. Comparison of the theoretical continuum and 
photometric measur ements (symbols) in different bands: B-, V-, R- and I-band photom etry from 



Massev et al. 



(120061 crosses); J-, H- and K-ba nd photometry f rom 



and Spitzer IRAC 3.6 and 4.5 /im fluxes from 



McOuinn et all (12007 



Cioni et al. 



(2009 



squares), 



, circles). The theoretical 



continua predicted by FASTWIND (solid line-best model, dash-dotte d line-enhanced M, see text 



for exp lanation) has been reddened following the extinction curye by 
Jl989 ). assuming Av/E(B-V) = 3.1. 



Cardelli. Clayton & Mathis 



-45- 




Fig. 11. — Color-magnitude and color-color diagrams constructed from the WFPC2/HST data, 
covering a field of 2.7x2.7 arcmin^ approximately centered at the location of NGC588. The 
reddening free color Q is defined as Q = {U —B) — 0.72 x (B — V). The two known Wolf-Rayet 
stars in NGC588 are identified by filled symbols, with the star marking UIT005. For reference, 
we also show an isochrone for a population ~5 Myrs old. 



